


Abstract.— Commercial and 
scientific bottom longline catches of 
alfonsino, Beryx splendens, from 
seamounts off New Caledonia were 
sampled to study length-frequency 
distributions. A total of 14,674 fish 
were measured. CPUE of Beryx 
splendens on two seamounts is mod- 
elled in terms of length and depth. 
The data show that mean length in- 
creases with depth; this is well de- 
scribed by a bivariate normal model 
that estimates catch for a given sea- 
mount. In addition, the data show 
that mean length also varies with 
the depth of the top of seamounts; 
this is described by a recursive 
model that is designed to predict ap- 
proximate catch for any seamount. 
The limitations of both models are 
discussed, particularly with regard 
to temporal variation. 
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A bottom longline fishery operated 
on the seamounts of the Exclusive 
Economic Zone (EEZ) of New Cale- 
donia from February 1988 to July 
1991.! Three vessels were involved 
but only one vessel was operated at 
any given time. The fishing effort, 
which totalled 4,691,635 hooks, fo- 
cused on five seamounts (B, C, D, 
J, and K) whose summits are lo- 
cated at depths ranging from 500 
to 750 m(Fig. 1). The target species, 
alfonsino, Beryx splendens, ac- 
counted for 92% of the catch by 
weight. This species has a world- 
wide distribution, from the equator 
to the temperate latitudes, and is 
fished by bottom trawl or longline. 
Alfonsino generally occupies waters 
between 200 and 800 m, although 
it has been caught at depths of only 
25 m and as deep as 1,240 m 
(Busakhin, 1982). Some authors 
have noted an increase in mean 
length with depth? (Yamamoto et 
al., 1978, Seki and Tagami, 1986), 
a trend which has been observed in 
other fishes (Heincke, 1913), par- 
ticularly some deep-water demersal 
species? + ? (Ralston and Williams, 
1988). There have been few studies 
relating the size distribution of 
alfonsino to depth. The objective of 
this paper is to describe an ap- 
proach for estimating the abun- 
dance of alfonsino by modelling its 
distribution in terms of fork length 
and depth of capture. À bivariate 
normal model describes this distri- 
bution for a given seamount and a 


recursive model predicts catch on 
any seamount. 


Material and methods 


Data 


Alfonsino were captured with long- 
line gear (Fig. 2). The main line, 
averaging 4,000 m, was held on the 
bottom by means of terminal an- 
chors and regularly spaced heavy 
sinkers that delimited five equal 
line sections. During a fishing trip 


! Grandperrin, R., and P. Lehodey. 1993. 
Etude de la pécherie de poissons profonds 
dans la zone économique de Nouvelle- 
Calédonie. Rapport final. Contrat de re- 
cherche ORSTOM / Territoire de Nouvelle- 
Calédonie. Nouméa: ORSTOM, Conv. Sci. 
Mer, Biol. Mar. 9, 321 p. 

Masuzawa, T., Y. Kurata, and K. Onishi. 
1975. Results of group study on population 
of demersal fishes in water from Sagami 
Bay to the southern Izu Islands— popula- 
tion ecology of Japanese alfonsin and other 
demersal fishes. Japan Aquatie Resources 
Conserv. Assoc. Fish. Res. Paper 28, 105 
p. [English translation held at Fisheries 
Research Centre Library, MAF, P.O Box 
297, Wellington]. 

Brouard, F., and R. Grandperrin. 1985. 
Deep bottom fishes of the outer reef slope 
in Vanuatu. South Pacific Commission 
17th Regional Technical Meeting on Fish- 
eries, W P 12, 127 p. 

Clark, M. R., and K. J. King. 1989. Deep- 
water fish resources off the North Island, 
New Zealand: results of a trawl survey, 
May 1985 to June 1986. N.Z. Fish. Tech. 
Rep. 11, 55 p. 

Dalzell, P., and G. L. Preston. 1992. Deep 
reef slope fishery resources of the South 
Pacific. South Pacific Comm. Inshore Fish. 
Res. Project. Tech. Doc. 2, 299 p. 


t 


ES 


in 


748 


Lehodey et al.: Modelling the distribution of Beryx splendens 


749 


* Tonga 


VT 
e E 
"om New Caledonia 


Australia 


New Zealand 


Depth (m) below 
surface 


0- 1000 
1000 - 2000 
2000 - 3000 

> 3000 


Figure 1 
Main seamounts fished for alfonsino, Beryx splendens, by bottom longline in New Caledonia. 





made by the longliner Humboldt from May to July 
1991 over seamounts B, C, D, J, and K,® the depth 
profile of the bottom was recorded on an echosounder 
as the line was set. The position and the depth at 
the exact time the terminal anchors and intermedi- 
ate sinkers were thrown overboard were also re- 
corded. The longliner Humboldt was equipped with 
a Doppler sonar current indicator which provided 
current velocity and direction at three selected 
depths. Data recorded suggest that the current ve- 
locity rapidly decreased with depth (Fig. 3A) and that 
horizontal drift was probably minimal.® On 23 occa- 
sions over the total of 73 longline sets, the depth of 
the bottom was recorded at the time the buoy was 
grabbed at the beginning of retrieval. This depth was 
compared with the depth of the corresponding ter- 
minal anchor recorded when the line was set. Depth 
difference was less than 10 m for 74 % of the paired 
comparisons (Fig. 3B) which indicates that either the 


ê Lehodey, P. 1991. Mission d'observations halieutiques sur le 
palangrier Humboldt. Campagne de péche du 30 mai au 12 
juillet 1991, Nouméa. ORSTOM Rapp. Missions, Sci. Mer Biol., 
Mar. 8, 44 p. 


drift of the line during sinking was limited or the 
slope of the bottom was slight. Therefore, despite the 
lack of a maximum depth recorder to determine the 
actual depth of the main line (Somerton and Kik- 
kawa, 1992), it was reasonable to assume that its 
configuration was similar to the depth profile indi- 
cated by the echosounder. 

The estimated depth of the sinkers was used to allo- 
cate a mean value of depth of capture d, = 1/2 (d, * d.) 
to all the fish caught on the same 800-m line section 
(Fig. 2). Ten meters, which is roughly half the length 
of the branch lines, was then added to each mean 
depth of capture d, to correct for bias introduced by 
the fact that catches may occur at any hook level. 
Figure 3C gives the depth variation within each sec- 
tion. Eighty-one percent of the variation in depth is 
less than 15 m and 92% is within the 0-25 m range. 
This indicates that in most cases the longline was 
nearly horizontal with the bottom. Therefore, the 
allocation of a single depth of capture to all fish 
caught on the same line section seems reasonable, 
particularly as the depth of capture data were ag- 
gregated into 25-m depth classes for analysis. Dur- 
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Figure 2 
Longline gear employed from the longliner Humboldt during a set (15) on 7 June 1991 
on seamount K (lat. 24'43'S; long. 170'06'E) and detail of a branch line. Main line is 
4,000 m long and divided into five sections (each section has 840 hooks [42 branch lines 
x 20 hooks] and is 800 m long). d, depth recorded on the echosounder at time t, when 
sinker was thrown overboard; d, = 1/2 (d, + d; , ,) + 10 = mean depth of capture allo- 
cated to all fish hooked on section 7. 








Figure 3 
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ing two other commercial cruises conducted by the 
longliners Hokko Maru and Fukuju Maru, observers 
recorded the maximum and minimum depths reached 
by the longline’: ? (Fig. 3D). 

Fork length (FL) was measured on a total of 14,674 
alfonsino. During the commercial fishing trips, fish 
to be measured were randomly sampled from each 
set. When there were only a few fish, they were all 
measured. As the samples varied in size between sets, 
the length-frequency distribution of all the alfonsino 
caught was estimated by multiplying the number in 
the sample by the ratio of the total number of 
alfonsino caught to the total number of alfonsino 
measured for each longline set. During several sci- 
entific cruises all of the fish were measured.! Because 


7 Anonymous. 1988. Rapport de la campagne de péche à la 
palangre profonde dans la zone économique de la Nouvelle- 
Calédonie. Hokko Maru 107: février-mai 1988. Territoire de 
Nouvelle-Calédonie, Service Territorial de la Marine Marchande 
et des Péches Maritimes, 57 p. 

8 Laboute, P. 1989. Mission d'observation halieutique sur le 
palangrier japonais Fukuju Maru du 21 nov. au 12 déc. 1988, 
Nouméa. ORSTOM Rapp. Missions, Sci. Mer Biol. Mar. 2, 15 p. 
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all the depth zones were not sampled in the same way, 
catch per unit of effort (CPUE=number of fish caught 
per 100,000 hooks) was taken as the abundance index. 

The data collected during the Humboldt cruise was 
used to model the distribution of CPUE in terms of 
length and depth over seamounts. The models were 
validated with data generated by the scientific 
cruises and with data collected on board the two ad- 
ditional commercial boats, Hokko Maru and Fukuju 
Maru. 'These commercial data are less precise be- 
cause only maximum and minimum longline depths 
were recorded. 


Preliminary data analysis 


Table 1 shows the mean fork length by depth zone 
for each seamount sampled during the Humboldt 
cruise. This table suggests a significant increase in 
mean length with depth on each seamount. In order 
to model this increase, the data must be as repre- 
sentative as possible of the fish population over its 
depth range. For this reason, only data recorded for 
seamounts B and J were used for modelling the in- 
crease in length with depth. 


Tabie 1 


Analysis of variance and multiple comparisons of mean fork length (cm) of alfonsino, Beryx splendens, sampled on 
five seamounts during the Humboldt fishing cruise (sample size in parentheses). 


Seamount 


Water depth of top B C 


(m) [502 m] [560 m] 


« 525 
525-549 
550-574 
575-599 
600-624 
625-649 
650-674 
675-699 
700-724 
725-749 
750-774 
775-799 

+ 800 

a* 0.0001 


33.27 (52) 
34.90 (665) 
37.11 (635) 
38.08 (82) 
40.39 (125) 


34.72 (405) 
36.78 (1083) 
38.95 (74) 


0.0001 





35.23 (262) 
36.24 (462) 
37.55 (913) 
38.16 (115) 


0.0001 





[630 m] 


0.0001 
0.0007 
0.0001 


35.73 (138) 
35.54 (295) 


0.0001 
0.0001 


40.60 (205) 


38.21 (309) 
39.90 (308) 
40.57 (53) 


0.4979 0.0001 





* [fa is less than 0.05 then the hypothesis that the means are the same in all classes is rejected. All individual means were compared pairwise with 
the multiple comparison test of Tukey-Kramer (in SAS, 1988). The shaded boxes indicate that the two included means are not significantly 


different at the 0.05 level. 
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CPUE distributions by size class and depth zone 
over seamounts B and J from the Humboldt are 
shown in Figure 4. A preliminary examination ofthe 
data revealed that they fitted portions of curves con- 
forming to a normal distribution. It was therefore 
assumed that for a given seamount the CPUE, in 
terms of length and depth, was distributed over a 
surface described by a bivariate normal distribution 
function delimited by the maximum and minimum 
of lengths and depths sampled. This assumption is 
the basis of the first modelling exercise (“bivariate 
normal model"). 

Tahle 1 also shows that, for a given absolute depth, 
mean length significantly decreases as the depth of 
the top of the seamount increases. This decrease sug- 
gests that the length distribution depends both on 
the absolute depth (in relation to the sea surface) 
and on the depth of the top of the seamount. Conse- 
quently, the bivariate normal model constructed for 
a given seamount may not be applicable to other sea- 
mounts whose summits lie at different depths. It is 
therefore necessary to construct a more general 
model (referred to as the "recursive model") which 
would predict extrapolated estimates of CPUE over 
any seamount hy taking into account both the abso- 
lute depth of the water column and the depth of the 
top of the seamount. Temporal validation of these 
two models requires data that were not used during 
model construction but were collected in the same 
area at different periods. Data collected on board RV 
Alis and the fishing vessels Hokko Maru and Fukuju 
Maru were used for model validation. 


Modelling method 


Bivariate normal model In the bivariate normal 
model, CPUE by length and depth are calculated on 
the basis of a bivariate normal distribution defined 
by the density function (1) 


1 
Bix, »Xg )z 


2120; Og lap 


| 1l l Xj — Hj | E (x, - Hi xa - Ha) (1) 
L O; 0; Og 


+ 





(Xa - Ha) } | 
653 P 


where x, is the length, x, is the depth, u, is the mean 
length, o, is the standard deviation of length, u4 is 
the mean depth, o, is the standard deviation of depth, 
andpis the regression coefficient of length on depth. 
Because sampling of the seamounts is limited up- 
wards by their summit (D.) and downwards by the 
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maximum depth accessible with the bottom longline 
(D, ), CPUE distributions will be modelled by a por- 
tion of the bivariate normal distribution (2) 


CPUE. (xj, Xg) =0 for x, END or eD. 
CPUE (xj, xq) = A Bix, x4) for D, S «D 


est a 


(2) 


where A represents the theoretical cumulative CPUE 
estimated over the field of definition of the entire 
bivariate normal distribution. The parameters À, 11, 
6, Hg 6, and p were estimated by a nonlinear regres- 
sion by using an iterative algorithm for sum of squares 
errors (SSE) minimization (SAS, 1988). 


Recursive model The recursive model should allow 
estimation of alfonsino CPUE by size class for sea- 
mounts for which no data are available except the 
depth of their summit. The size structure variation 
shown in Table 1 should be taken into account in the 
development of this model, i.e. 1) for a given sea- 
mount, mean length increases with depth and 2) for 
a given depth zone, mean length decreases as the 
depth of the top of the seamount increases. In theory, 
the distribution of a population of alfonsino over any 
seamount could then be taken as the superposition 
ofthe distributions of two subpopulations: one popu- 
lation would be influenced only by the absolute depth 
while the other would be influenced by the depth of 
the top of the seamount. This model attempts to ex- 
plain how the fish population of a given seamount 
would in theory redistribute itself if it were to mi- 
grate and settle on another seamount. Consider a 
hypothetical population whose cumulative CPUE 0.) 
by length and depth is distributed over its seamount 
of origin according to a bivariate normal distribu- 
tion function of unknown 41, O, Hg Og, and p param- 
eters. For the estimation of these parameters, the 
top of the hypothetical original seamount will be as- 
sumed to be exactly level with the sea surface in or- 
der to include the entire depth zone that could be 
inhabited by alfonsino. The new CPUE distribution 
on seamounts jJ, j*1, etc. . . (Fig. 5) will depend on 
the initial parameters of the distribution over the 
original seamount as well as on parameter p, the 
probability that the fish will redistribute according 
to absolute depth, and ( 1-p), the probability that the 
fish will redistribute according to the depth of the 
top of the seamount. At each “leap” to a deeper sea- 
mount (from j to j4 1), the subpopulation inhabiting 
a given depth zone D, will split into two groups: one 
group will stay in the same depth zone D, with a prob- 
ability p (or will migrate elsewhere if this zone is no 
longer available on the new seamount) and the other 
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Figure 4 


Actual CPUE for alfonsino, Beryx splendens, by fork length (cm) and depth (m), recorded on seamounts B 
and J during the fishing cruise carried out by the longliner Humboldt, and predicted CPUE from the 
bivariate normal model and the recursive model. 
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Figure 5 
Changes in distribution of alfonsino, Beryx splendens, according to the 
recursive model of two subpopulations Y and Z that migrate from a sea- 
mount j to deeper seamounts j+1, j+2, /+3. Y=subpopulation distributed 
on a seamount j at absolute depth D; Z=subpopulation distributed on a 
seamount j at absolute depth D, ,; p-probability that fish will redistrib- 
ute according to absolute depth; 1-p=probahility that fish will redistrib- 
ute according to the depth of the top of the seamount; *=original sea- 


mount. 





group will move down to zone D,,, with a probability 
l-p (Fig. 5). Thus, it is possible to determine the 
CPUE (X, ,,) for the population of a depth zone D, 
on a seamount /, for a size class k in terms of the 
subpopulations of zones D, and D, , on the higher- 
level seamount j-1. This is expressed as follows: 
X, = PX, 


Lit a 


pb PAS (3) 


Specifically, if p=0, CPUE is distributed solely ac- 
cording to the depth of the top of the seamount and 
if p=1, it is distributed solely according to the abso- 
lute depth. If the parameters Àp Hp Op Ly, Og p. and 
p. estimated from a known seamount length-depth 
distribution are known, it is possible to calculate all 
the CPUE (X, p?) values for any seamount j (deeper 
or shallower), depth zone D, and size class k. The 
foregoing seven parameters can be estimated by mini- 
mizing the SSE between the CPUE recorded on one 
of the best sampled seamounts (B or J) and the CPUE 
estimated by Equation 3. This estimation 1s per- 
formed by a nonlinear regression (SAS, 1988). 


Results 


Bivariate normal mode! 


Application of a bivariate normal model implies that 
mean length can be deduced from depth by a linear 
regression weighted by the CPUE x, 2 ax, + b where 
a and b are constants). The results of this regression 
for seamounts B and J show mean length and depth 
to be significantly correlated ( Table 2). Consequently, 
the bivariate normal model can be tested for each of 
these seamounts. 

The parameters of the bivariate normal model were 
calculated separately for seamounts B and J ( Table 
3). The determination coefficient,? R?, for seamounts 
B and J respectively equals 0.87 and 0.93. The re- 
sidual analysis was carried out to test the fit of the 
model to the data from the Humboldt cruise on sea- 
mounts B and J. The results show the residuals are 


2 


SR s [xu - vu vl with Y = CPUE. 
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Table 2 
Bivariate normal model: CPUE - weighted linear regression of length of alfonsino, Beryx splendens, on depth. 


No. of fish Min. depth Max. depth 
Seamount measured (m) (m) P a Hy: 7 =0 b 


B 1,557 516 615 0.549 < 0.0001 rejected 0.037 
J 1,957 606 761 0.486 < 0.0001 rejected 13.122 





p = regression coefficient. 

a = significance probability of the regression under the null hypothesis that the statistic is zero. 
H, = null hypothesis i.e. length and depth are independent, if a < 0.05, H, is rejected. 

a and b = parameters of the linear regression. 





satisfactory; in particular, the residuals are centered 
on zero, are not correlated with the length and depth 
variables, and have a constant variance (Table 4; Fig. 
6). These characteristics indicate a good fit of the 


Table 3 


Bivariate normal model: predicted parameters for 
seamounts B and J. SD=Standard deviation. 


bivariate normal model to the data as demonstrated eM eat Eee 
by comparison of actual and predicted CPUE (Fig. 4). 
Extrapolation of the model to the data not used in B J 


the modelling exercise is unsatisfactory because the 1,557 1,957 


mean value of the residuals is not centered on zero 


d . Parameters Estimation SD Estimation SD 
for the Fukuju Maru data and because the residuals 


are correlated with the length variable for the RV À 3.66x109  0.85x109 0.5x109° 4.5x10° 
Alis data and with the depth variable for the Hokko B 35.22 0.96 12.0 52.0 
Maru data (Table 4; Fig. 6). This suggests the exist- Ha 542.87 14.49 -123.1 1,666.7 
ence of factors affecting the population’s distribution = 591 191 d a6 

o, 71.07 17.04 272.5 317.1 
not accounted for by the model. p 0.75 0.10 0.9 0.2 





A-theoretical cumulative CPUE. 
A u,2mean length (cm). 
Recursive model ji-mean depth (m). 
o=standard deviation of length. 
: 3 o,=standard deviation of depth. 
The parameters of the recursive model were esti- p?-regression coefficient of length on depth. 


mated separately for seamounts B and J. The deter- 





Table 4 


Bivariate normal model: results of analysis of residuals (£) for fit control and temporal validation of the model for 
seamounts B and J. 








No. of fish 
Cruise Seamount measured He =0 B Hosp Oy 


Fit control Humboldt 1,557 not rejected not rejected 0.149 not rejected 
Humboldt 1,957 not rejected not rejected 0.431 . not rejected 


Temporal 
validation Hokko Maru 
RV Alis 


2,840 0.262 not rejected 0.601 notrejected | «0.0001 rejected 
1,688 0.908 not rejected 0.016 rejected 0.391 not rejected 


B 
B 
Fukuju Maru J 4,320 0.0002 rejected 0.265 notrejected 0.284 not rejected 








H, : € = 0. The mean value of the deviations between estimated and observed CPUE is 0, If a, is «0.05, H, is rejected 
P,7 regression coefficient of € on length. 

p,-regression coefficient of € on depth. 

Hyp, = 0. If o, is« 5%, Hy is rejected. 

Hy; = 0. If a, is < 5%, H, is rejected. 
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Figure 6 


Bivariate normal model and recursive model for seamounts B and J: distribution of predicted CPUE of alfonsino, 
Beryx splendens, in relation to actual CPUE and distributions of residuals in relation to length (cm) and depth ( m). 


Dotted lines delimit the confidence interval at a=0.05. 





Table 5 


Recursive model: estimated parameters of the dis- 
tribution of CPUE on the hypothetical original sea- 
mount calculated from seamount J data. SD = stan- 
dard deviation. 





Parameters Estimation SD 





22.86 91.3 
-66.17 499.9 
5.98 16.5 
30.27 9917 
0.87 0.70 
0.09 0.04 
1.15x10^ 2.33x10^ 








u,2mean length (cm: 
1,2mean depth im! 
6,=standard deviation of length 
,zstandard deviation of depth 
egression coefficient of length on depth 
p=probability that fish will redistribute according to absolute depth 
4 ztheoretical cumulative CPUE 








mination coefficient, R?, calculated for seamount J 
equals 0.82, while for seamount B it equals 0.69. 
Therefore, the parameters estimated for seamount 
J were used in the model (Table 5). The residuals 
resulting from fitting the model to data from the 
Humboldt cruise on seamount J are satisfactory; in 
particular, they are centered on zero and are not cor- 
related with the studied variables (Table 6; Fig. 6). 
These features indicate a good fit of the recursive 
model to the data as demonstrated by comparison of 
actual and predicted CPUE for seamount J (Fig. 4). 
It is interesting to note the low value of p (close to 
0.1 as shown in Table 5), which indicates that the 
seamount top depth parameter has greater impact 
on the length distribution than does the absolute 
depth parameter. 

Spatial validation was carried out for seamount B 
from the data collected during the Humboldt cruise 
( Table 6). The residuals are centered on zero and not 
correlated with the length and depth variables. How- 
ever, since their variance is not constant ( Fig. 6) and 
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Table 6 
Recursive model: results of analysis of residuals (£) for fit control and validation of the model for seamounts B and J. 


No. of fish 


Cruise Seamount measured wi Hoe = 0 H,:p,=9 Oy H,=p,=0 


Fit control Humboldt J 1.957 0.87 not rejected 0.85 not rejected 0.06 not rejected 


Spatial 
validation Humboldt 1-557 0.13 not rejected 0.06 not rejected 0.99 not rejected 
Temporal 
validation Hokko Maru B 2,840 0.15 
RV Alis B 1,688 0.01 
J 


Fukuju Maru 4,301 0.007 


not rejected 0.76 
rejected 0.005 
rejected 0.67 


not rejected <0.0001 
rejected 0.35 
not rejected 0.08 


rejected 
not rejected 
not rejected 


Hy : € = 0. The mean value of the residuals is 0. If a, is <5“, H, is rejected. 
p, = regression coefficient of £ on length. 
p; = regression coefficient of € on depth. 
H :p, = 0. If a, is <5%, H, is rejected. 
Hp, = 0. Ifa, is <5%, Hy is rejected 


since the standard deviations of the parameters are 
high (Table 5), spatial extrapolation of the model to 
seamount B is rather crude as demonstrated by com- 
parison of actual and predicted CPUE (Fig. 4). 

Temporal validation was carried out on data from 
the Fukuju Maru and RV Alis fishing on seamounts 
B and J (Table 6). It is unsatisfactory because the 
mean values ofthe residuals are not centered on zero 
and the residuals are correlated with the length vari- 
able for the RV Alis data and with the depth vari- 
able for the Hokko Maru data (Fig. 6). As with the 
bivariate normal model, this suggests the existence 
of factors not accounted for by the model. 


Discussion 


Alfonsino length structure variation observed over 
the seamounts of New Caledonia is similar to that 
noted in Japan? and in New Zealand" (Massey and 
Horn, 1990) where it was assigned to age-specific 
migrations. In Japan, it was noted that alfonsino 
move south as they grow,” young fish predominate 
over some seamounts and old fish predominate over 
other seamounts. In New Caledonia, age segregation 
over the seamounts is so marked that it has been 
possible to describe it mathematically. 

The bivariate normal and recursive models appear 
to be complementary. The bivariate normal model 
provides an instantaneous picture of alfonsino popu- 
lation distribution on a given seamount; it provides 
good CPUE estimates provided a sufficient amount 


1? Horn, P. L., and B. R. Massey. 1989. Biology and abundance of 
alfonsino and bluenose off the lower east coast, North Island, 
New Zealand. N. Z. Fish. Tech. Rep. 15, 31 p. 





of length and depth data are available. The recur- 
sive model takes into account the dynamic nature of 
the population's distribution as it allows the extrapo- 
lation of CPUE obtained for one seamount to sea- 
mounts that were not sampled. It allows preliminary 
population estimation of unexploited stocks. Depend- 
ing on current economic parameters, the model might 
be used to indicate the depths at which fishing is 
most economic. Once a fishery is operational, more 
refined data will be available, which will enable the 
bivariate normal mode! to be applied and stock man- 
agement parameters defined for each of the sea- 
mounts fished. 

The poor results obtained for the temporal valida- 
tion could be due to poor precision of the depth data 
collected from the longliners Hokko Maru and Fukuju 
Maru. Also, neither of the models incorporate a time 
factor. The data were collected from cruises carried 
out in different years and in different seasons. Hence, 
it is unlikely that conditions remained stable, par- 
ticularly with regard to exploitation history, repro- 
ductive behavior, or long-term climatic variations. 

Fishing methods and strategies were not modified 
during the fishing period considered. Therefore, the 
catches are probably representative of the standing 
stock of alfonsino within the size limits determined 
by the selectivity of the fishing gear. Since the daily 
observation window did not change, vertical trophic 
migrations would seem unlikely to contribute to the 
observed variability. With regard to sex as a source 
of variability, although the mean length of females 
exceeds that of males!! (Kotlyar, 1987; Massey and 


11 Lehodey, P. 1994. Les monts sous-marins de Nouvelle-Calédonie 
et leurs ressources halieutiques. Thése de doctorat de 
l'Université Francaise du Pacifique, 398 p. 
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Horn, 1990), Humphreys et al. (1984) have shown 
that sexual dimorphism is not responsible for the 
existence of different size groups of alfonsino. Marked 
declines in CPUE are observed in the Southern Hemi- 
sphere during summer. This season corresponds to 
the alfonsino breeding period in New Caledonian 
waters.!! The summer decline in catch rate could be 
due to breeding migrations drawing the fish to 
spawning grounds that are different from the fish- 
ing grounds? (Chikuni, 1971) or to changes in vul- 
nerability to the gear owing to seasonal physiologi- 
cal or behavioral changes (Ricker, 1980). Data used 
to build the models were collected on board the 
Humboldt during the winter season. Data used to 
validate the models were collected on board Fukuju 
Maru and Hokko Maru at the beginning and end of 
the warm season and during six scientific cruises, 
five of which were carried out in summer. This sug- 
gests that reproductive seasonality might be a fac- 
tor in the poor temporal validation of the models. 
Other sources of temporal variation might be re- 
lated to the environment. The ocean habitat of 
alfonsino is not affected by continental influences but 
is subject to hydrological fluctuations affecting the 
deep-water masses. Some of these influences are of 
short period such as internal waves and tidal cur- 
rents (Eriksen, 1985; Roden, 1987), whereas others 
recur at longer intervals such as seasonal variations 
in ocean currents and multi-annual hydroclimatic 
anomalies of the El Nino Southern Oscillation 
(ENSO) (Delcroix and Hénin, 1989). Such fluctua- 
tions might have an impact on alfonsino stock struc- 
ture, either at the recruitment stage (survival and 
dispersal of eggs and larvae) or by modification of 
the behavior of adults (migrations from one seamount 
to another). However, it is difficult to demonstrate 
the effect of these fluctuations on the presence and 
catchability of fish. It is even more difficult to ex- 
plain the very large differences in fishery productiv- 
ity observed between seamounts of identical depth, 
located only a few dozen miles apart and appearing 
to have the same hydrological environment. Seafloor 
topography and bottom type might account for these 
differences, but other hypotheses can be postulated, 
some based on the existence of a low-energy hydro- 
thermalism (Rougerie and Wauthy, 1990) and oth- 
ers on a hydrological anomaly called "Taylor's col- 
umn," which could enhance species sedentarity 
(Royer, 1978; Genin and Boehlert, 1985; Roden, 1987; 
Dower et al., 1992; Sime-Ngando et al., 1992). Fluc- 
tuations in intensity of this anomaly, or its disap- 
pearance, could also be responsible for the variations 
in productivity observed over time over a given sea- 
mount (Boehlert and Genin, 1987). These unknown 
environmental fluctuations cause problems in the 


interpretation of results from exploratory and com- 
mercial fishing cruises carried out over seamounts. 
The data collected at a given location constitute an 
instant picture of a stock whose abundance is likely 
to vary, irrespective of fishing effort, as a result of 
unknown environmental variations. In other words, 
the fertility of the seamounts could vary quite un- 
predictably over the history of a fishery. Conse- 
quently, modelling the distribution of a stock should 
be confined to a relatively small temporal sampling 
scale. 


Conclusion 


The bivariate normal model and the recursive model 
provide complementary interpretations of length dis- 
tribution in terms of depth of alfonsino fished on the 
seamounts of New Caledonia by the bottom longline 
fishery. They could be useful for the proper manage- 
ment of fisheries over seamounts, where stocks are 
known to be vulnerable (Sasaki, 1986) because of the 
limited habitat afforded by seamounts and the slow 
growth rate of deep-water species. However, it would 
appear that annual or seasonal factors, in particu- 
lar those which account for recruitment fluctuations 
and behavioral changes linked to reproduction, will 
need to be incorporated into the models before they 
can be generalized. A better understanding of the 
functioning of the ecosystems concerned would also 
assist in establishing the limits of generalization, 
particularly with regard to depth and area inhab- 
ited by alfonsino. These models could possibly be 
adapted to other deep-water species such as certain 
snappers and groupers. 
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